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Transport modeling of idealized, cone-guided fast ignition targets indicates the severe challenge posed by 
fast-electron source divergence. The hybrid particle-in-cell [PIC] code Zuma is run in tandem with the 
radiation-hydrodynamics code Hydra to model fast-electron propagation, fuel heating, and thermonuclear 
burn. The fast electron source is based on a 3D explicit-PIC laser-plasma simulation with the PSC code. 
This shows a quasi two-temperature energy spectrum, and a divergent angle spectrum (average velocity- 
space polar angle of 52°). Transport simulations with the PIC-based divergence do not ignite for > 1 MJ 
of fast-electron energy, for a modest (70 fim) standoff distance from fast-electron injection to the dense fuel. 
However, artificially collimating the source gives an ignition energy of 132 kJ. To mitigate the divergence, 
we consider imposed axial magnetic fields. Uniform fields ~50 MG are sufficient to recover the artificially 
collimated ignition energy. Experiments at the Omega laser facility have generated fields of this magnitude 
by imploding a capsule in seed fields of 50-100 kG. Such imploded fields are however more compressed in 
the transport region than in the laser absorption region. When fast electrons encounter increasing field 
strength, magnetic mirroring can reflect a substantial fraction of them and reduce coupling to the fuel. A 
hollow magnetic pipe, which peaks at a finite radius, is presented as one field configuration which circumvents 
mirroring. 

PACS numbers: 52.57.Kk, 52.65.Kj, 52.65.Rr, 52.65.Yy 



I. INTRODUCTION 

The fast ignition approach to inertial fusion exploits a 
short-pulse, ultra-intense laser to heat an isochoric hot 
spot to ignition conditions 1 . Unlike the central hot- 
spot approach, fast ignition separates dense fuel assem- 
bly from hot-spot formation 2 . This opens the prospect 
of high energy gain with less laser energy, and may be 
an attractive avenue for inertial fusion energy. The first 
integrated but sub-ignition scale experiments were per- 
formed at Osaka in 2001-2002, and explored the cone-in- 
shell geometry 34 . Subsequent similar experiments were 
done at VulcarP, Omega EP^, and Osaka in 2009 and 
2Q2(jZHl The 2002 Osaka experiments were interpreted 
to show high coupling of short-pulse laser energy to the 
fuel, of order 20%. All the later experiments show lower 
coupling, the best being 10-20% coupling in the 2010 
Osaka experiments^. This work suggests the reduced 
coupling seen in 2009 at Osaka was due to higher pre- 
pulse in the short-pulse laser. Pre-pulse energy creates 
an underdense pre-plasma in which the laser converts to 
over-energetic electrons, and this source is farther away 
from the fuel. The negative impact of pre-plasma on fast- 
electron generation inside a cone has been reported, e.g., 
in Refs. I9HTTI Coupling efficiences at small scale do not 
directly apply at ignition scale. 

This paper presents integrated fast-ignition modeling 
studies at ignition scale, which is well beyond parameters 
currently accessible by experiment. We utilize a new, hy- 
brid PIC code ZumsP^ to model fast electron transport 
through a collisional plasma, with self-consistent return 
current and electric and magnetic field generation. To 
alleviate the need to resolve light waves or background 



Langmuir waves, Zuma does not include the displacement 
current in Ampere's law, and employs an Ohm's law (ob- 
tained from the inertialess limit of the background elec- 
tron momentum equation) to find the electric field. We 
recently coupled Zuma to the radiation-hydrodynamics 
code HydrsP^, which has been widely used to model in- 
ertial fusion and other high-energy-density systems. 

We do not model the short-pulse laser, but instead 
inject electrons with a specified distribution into Zuma. 
The source electron spectrum is a key element of this 
approach. We o btain the spectrum by using the particle- 
in-cell code PSCP^ to perform a 3D full-PIC simula- 
tion of the laser-plasma interaction (LPI). This gives a 
quasi two-temperature energy spectrum, with a (cold, 
hot) temperature of (19, 130)% of the so-called pondero- 
motive temperature as defined below^ at the nominal 
laser intensity. It is generally seen in PIC simulations 
that LPI occurring at lower density produces more ener- 
getic electrons. The experimental understanding of fast 
electron energy spectra is not entirely clear. Ma et al. re- 
cently reported experimental evidence indicating a two- 
temperature energy spectrum^, albeit at lower energies 
and shorter pulses than considered here for ignition. 

The PIC-based angle spectrum is very divergent, with 
an average polar angle in velocity space of 52° or an inte- 
grated solid angle of 4.85 sterad. A large divergence has 
been reported in other PIC simulations, such as Refs. IT71 
and [18| Experimental evidence for a significant diver- 
gence comes from modeling by Honrubia et al. 19 of K a 
data obtained by Stephens et alW^. More recent work 
by Westover et al. 21 also indicates a substantial source 
divergence. 

Our Zuma-Hydra modeling with a realistic fast elec- 
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tron source (both energy and angle spectra) and an ide- 
alized fuel mass located 70 /im from the electron source 
indicates poor coupling to the fuel hot spot, with > 1 MJ 
of fast electrons inadequate to ignite. Artificially col- 
limating the electron source dramatically improves the 
picture, with an ignition energy of 132 kJ. This is much 
higher than the ideal estimate of 8.7 kJ absorbed in the 
ignition hot spot (detailed below), due largely to the en- 
ergy spectrum being too hot for the electrons to stop fully 
in the hot spot. In Ref. [22j we report in more detail the 
effects of the energy spectrum, as well as E and B fields, 
on the ignition requirements for an artificially-collimated 
fast electron source. We merely note here that, for our 
particular plasma condition profiles, using the complete 
Ohm's law Eq. (14) reduces the fast-electron coupling 



to the fuel compared to the case of no E or B fields, 
while using the resistive Ohm's law E = 77J5 increases 
the coupling over the no-field case. This is likely due 
to Vn e x VT e magnetic fields that develop at the outer 
radius of the dense fuel and push the fast electrons to 
larger radius, as observed earlier in Ref. [23j 

The focus of this paper is on mitigating the beam di- 
vergence by imposed magnetic fields. In particular, we 
do not pursue here other attractive options, such as field 
generati on by resistivity gradients 2 ^. Cylindrical and 
spherica l 26 ^ 27 ! implosions at the Omega laser have com- 
pressed seed magnetic fields ~10 kG to strengths of 20-40 
MG. We show that a uniform, initial axial field of 50 MG 
almost recovers the ignition energy of the artificially col- 
limated beam. We stress that such magnetic fields do 
not collimate the electrons (that is, reduce the velocity- 
space divergence), but rather confine them in space: once 
they emerge from the confining field, they still have their 
initial divergence. A magnetic field that increases in the 
axial direction leads however to substantial reflection due 
to magnetic mirroring. We explore a hollow magnetic 
pipe, that peaks at a finite radius, as one method to cir- 
cumvent mirroring. Pipes with a peak field of 50 MG and 
radial FWHM of (20, 30) /im ignite for (211, 158) kJ of 
fast electrons, compared with 132 kJ for an artificially- 
collimated source with no imposed field. 

The paper is organized as follows. Sec. II describes the 
fast electron source derived from a full PIC simulation of 
the LPI. In Sec. Ill we detail the Zuma model, and how 
it is coupled to Hydra. Sec. IV presents Zuma- Hydra 
results with realistic fast electron divergence and an arti- 
ficially collimated source. We study ways to mitigate the 
source divergence with several imposed magnetic fields 
(uniform, increasing, and a hollow pipe) in Sec. V. In 
Sec. VI we make some concluding remarks. 



II. ELECTRON SOURCE FROM FULL-PIC 
LASER-PLASMA MODELING 

The distribution of fast electrons produced by the 
short-pulse LPI is a crucial element of fast ignition. The 
computational scale of integrated LPI-transport-hydro 



simulations is currently prohibitive. We excite a fast 
electron source derived from full-PIC LPI simulations in 
our transport modeling. This neglects feedback of the 
transport on the LPI, e.g. the exact details of the return 
current, and also loses some detail in the fast-electron 
source. It does call out the major design challenges, and 
allows for the development of ideas to mitigate them. 

Our LPI simulations are performed with the relativistic 
PIC code PSCP^. This code has recently been extended 
to include a hybrid model, valid for collisional plasmas^. 
The results presented here do not use the hybrid model, 
and are all explicit, full-PIC calculations, with the com- 
plete Maxwell equations and no collision operator. We 
note that more recent PSC simulations of several ps dura- 
tion show the convergence of profiles with and without an 
initial pre-plasma, and the development of a third, super- 
hot component (temperature several times the pondero- 
motive value given below) besides the two components 
present in the results discussed here. Our source will be 
updated with these new results in the near future. 

The specific PSC run used for the fast electron source 
was as follows. The geometry was 3D Cartesian, and 
the electron density at time 360 fs over a 2D plane is 
shown in Fig. [I] The domain extended from —30 to 
+30 /im in the two transverse (x and y) directions and 
from to 40 /im in z (nominal direction of laser propaga- 
tion). The particle and field boundary conditions (BCs) 
were periodic in x and y. In z, the particle BCs were 
thermalizing re-emission, while the field BCs were ra- 
diative (outgoing- wave). The initial plasma profile was 
n e /n cr = 100 for z > zq and exp[(z — zo)/3.5/im] for 
z < zo with zo =20 /im and n cr = 1.115 • 10 21 /Ao[/im] 
cm -3 , the non-relativistic laser critical density. This pro- 
file was chosen to replicate the pre-plasma produced by a 
small pre-pulse (~1-10 mJ) in the short-pulse laser (e.g. 
growing from ASE). Both electrons and deuterium ions 
(Z/A = 1/2) were mobile. The uniform cell size was 
Ax = Ay = Ao/16 and Az = Ao/16.375. The time step 
was cAt = 0.421 Ax. There were (twelve, four) numeri- 
cal macro-particles per cell for (electrons, ions). The run 
required about 160k cpu-hours to complete. 

The laser had a vacuum wavelength of Ao =1 /im and 
a vacuum focal spot at z =10 /im with radial intensity 
profile I(r) = Iql exp [— (r/18.3 /im) 8 ] . For a given laser 
power and maximum (hard-edge) spot radius, a flat as 
opposed to peaked (e.g., Gaussian) profile reduces the 
average intensity and gives a cooler spectrum. We chose 
Iql = 1.37 • 10 20 W/cm 2 , corresponding to a normalized 
vector potential ao = 10. The ponderomotive tempera- 
ture, as defined in Ref. [16j for the peak laser intensity is 

T p /m e c 2 = [1 + al] 1/2 - 1 = 9.05 , or T p = 4.62 MeV. 
We simply use T p to denote an energy scale, without any 
implication about what the fast electron distribution is 
(which is discussed below). The laser pulse was ramped 
up to peak intensity over 30 fs. 

We took all electrons at the time 360 fs, in a cylindri- 
cal "extraction box" from z = 20 — 25 /im and radius 
30 /im. This box is deep enough into the overdense re- 



°n i i r o 

20 ,^ 40 60 

y A 

FIG. 1. (Color online) Electron density at time 360 fs in the 
PSC run used to characterize the fast electron source. The 
white box indicates the extraction box, and the black box 
indicates the source box in the hybrid-implicit LSP run. The 
laser was incident from z — with a vacuum focus at z — 10 
/im. 



gion that the laser did not propagate there. We also 
selected only electrons with 0.55 < E [MeV] < 29.5 (or 
0.12 < E/T p < 6.37) and v z > 0, where E = m e c 2 (j-l). 
This is done to eliminate the return current and back- 
ground heating. Some of this heating is unphysical grid 
heating, and some is a legitimate kinetic energy transfer 
between the fast electrons and the background plasma at 
the 100n cr density used in the PIC simulation for numer- 
ical reasons. This heating is expected to play a negligible 
role at the densities assumed in our hybrid simulation^. 

For our transport studies, we do not inject a fast elec- 
tron source in an analogous extraction box. In partic- 
ular, such a source would need a radially outward drift 
that varies with radius. Instead, we excite fast electrons 
in a "source box" analogous to the laser absorption re- 
gion, such that after propagating a small distance into 
an equivalent extraction box, the transport-code electron 
distribution matches that of the PIC electrons. This 
method automatically handles a host of issues regard- 
ing propagation from the source to extraction regions 
(e.g., finite "view factors" that vary with angle and ra- 
dius), and provides the overall laser to electron conver- 
sion efficiency. The source electron intensity^! is acEl(r) 
where I(r) is the vacuum laser intensity given above, and 
&ce — 0.52 is an overall laser-to-electron power conver- 
sion efficiency, olce was chosen so the total fast elec- 
tron kinetic energy in the PSC and transport-code ex- 
traction boxes match. The source intensity is varied in 
space and time only by varying the rate of excitation, not 
the velocity-space distribution. 

To arrive at the transport-code distribution excited in 
the source box, we performed a hybrid implicit-PIC sim- 
ulation with the LSP plasma simulation code^, with ki- 
netic fast electrons and fluid background species (this is 
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FIG. 2. Source energy spectrum E fs (solid) and its running 
integral (dash) for the analytic form in Eq. |l]). The circles 
indicate the limits of the domain taken in the extraction box, 
and injected into our Zuma-Hydra simulations. 

the only time LSP is used in this paper). The LSP source 
box was located from z = 10 to 15 /im, and the plasma 
was uniform 10 g/cm 3 carbon at 100 eV. We used this 
much denser background than the PSC simulation be- 
cause the fluid model is only valid at high collisionality, 
and our transport studies are performed in compressed 
matter. The difference between free-particle propagation 
and the full LSP results are small, indicating that forces 
are not important as electrons transit from source to ex- 
traction boxes. 

The following LSP source gave electrons in the extrac- 
tion box that agreed adequately with the PSC extrac- 
tion box electrons. The source velocity-space distribu- 
tion is azimuthally symmetric and given by d 2 N/dEdO = 
NofE(E)fo(6), so that (fE,fe) are proportional to the 
ID distributions (dN / dE , dN / dO) . No is an overall nor- 
malization factor. tan# = [v z / (v 2 + v 2 ) x / 2 ] defines the 
polar angle in velocity space. For PSC runs at ignition 
powers and wide focal spots, we generally find the angle 
spectrum does not vary much with energy (see Fig. [3|. 
This justifies our ID factorization; the method can be 
easily extended to several energy bins each with different 
fe- ' 

A. Fast electron energy spectrum 

For /e we use the ID energy spectrum dN/dE found 
in the PSC extraction box. This is well-fit by a quasi 
two-temperature form: 

fs(e) = - exp [-e/n] + 0.82 exp [-e/r 2 ] . (1) 

e = E/T p is the ponderomotively scaled energy, and we 
assume as we vary the laser intensity and wavelength that 
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FIG. 3. (Color online) Average velocity- space polar angle 
vs. electron kinetic energy in the extraction box for the PSC 
full-PIC (solid black) and LSP implicit-PIC (solid red) runs. 
The classical ejection angle O c from Eq. (pi is plotted in 
dashed blue. 



dN/dE scales in this manner. Figure [2] plots this analytic 
form, as well as its running integral. The temperature- 
like parameters have the values t\ = 0.19 and T2 = 1.3. 
These correspond to a relatively "cold" component pro- 
duced by LPI near or above n cr , and a slightly hotter 
than T p component arising from underdense LPI. The 
1/e factor on the cold part improves the fit at low energy, 
although this may change as better ways to eliminate re- 
turn current and background heating are developed. We 
only inject over the domain 0.12 < e < 6.37 , which is the 
domain taken from the PSC extraction box. The average 
injected electron energy is (e) = 1.02, while only 24% of 
the injected energy is carried by electrons with e < 1. 
This is unfortunate for ignition with the laser intensities 



we contemplate, as discussed in Sec. IV 



B. Fast electron angle spectrum 

The average angle in the extraction box as a function 
of electron energy is displayed in Fig. [3J The PSC full- 
PIC becomes slightly more collimated at higher ener- 
gies, while the LSP implicit-PIC is essentially indepen- 
dent of energy. The agreement is excellent for E < T p , 
but the LSP 6 is slightly larger at large energies. We con- 
sider the energy dependence of to be weak enough to 
ignore and use a factorized source. Both PIC simulations 
have much larger and much less decrease with energy, 
than the classical ejection angle C given by^ 
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FIG. 4. (Color online) Solid angle spectra in extraction box 
for PSC run as described in text (black), LSP run extraction 
box (red), and source fa from Eq. Q with an arbitrary 
scale factor (blue dashed). The source spectrum is broader 
than those in the extraction boxes due to limited view factor 
at large angles. Note that the first two curves are given in 
physical units. 



This result obtains for a single electron in a focused laser 
beam in vacuum, not including plasma effects. 
The source angle spectrum we use is 



fo(0) = 2<irsm(0)fa fa = exp - (0/AO) 



(3) 



Q represents solid angle, related to by dQ = 2tt sm(9)d9. 
The value of the parameter AO that gives good agree- 
ment with the angle spectrum in the extraction box is 
AO = 90° . Figure [4] displays the source fa as well as 
the angle spectra in the PSC and LSP extraction boxes. 
The resulting extraction box angle spectrum is somewhat 
narrower than this source, due to limited view factor at 
large 0. In addition, the LSP extraction spectrum is de- 
pleted at small angle compared to the PSC spectrum, and 
may slightly overstate the divergence (although there are 
few particles at small due to the sin# Jacobian). We 
stress that AO is a parameter in a function, and not an 
observable quantity. The average which has physical 



meaning, is 



Io /2 dOf90 ^ T[3/4] 



(0) 



Io /2 defo 



r l/2 



A0 = 0.691A6> A0 < 1. 



(4) 



T is the mathematical gamma function. (0) and the rms 
are shown vs. A# in Fig. [5j Note that for large A0, 
(0) falls below the approximate linear result given above. 
For AO = 90° we find (9) = 52° and rms = 56°. The 
integrated solid angle Qi is J Q 27r dVt fa = 4.85 sterad (for 
an isotropic distribution, fa = 1 and = 2tt). This is 
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a substantial divergence; the rest of this paper is focused 
on its consequences for fast ignition requirements, and 
mitigation ideas based on imposed magnetic fields. 

We briefly note that our source is symmetric in az- 
imuthal angle. However, it quickly develops a radially 
outward drift as it propagates. That is, the average an- 
gle, between an electron's position and velocity vec- 
tors decreases. For free motion, in the far-field limit they 
become parallel: r ~ vt. In the extraction box the LSP 
and PSC electrons have similar x distributions vs. radius. 
Debayle^ has recently discussed the role of an intrinsic 
radial drift produced by a Gaussian laser; we find such a 
drift naturally develops due to propagation from a sym- 
metric source. 
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electric field is found from Ohm's law (i.e., the momen- 
tum equation for the background electrons in the limit 
m e dv e b/dt — >• 0), and the background return current is 
found from Ampere's law without displacement current. 
This reduced-model approach is similar to GremilletP^, 
HonrubisP, DavieeP, and Cohen et al. 28 (although the 
last approach uses particles to describe the collisional, 
"fluid" background). This combination eliminates both 
light and Langmuir waves, and allows stable modeling 
of dense plasmas without needing to resolve these fast 
modes. An alternative approach to dense-plasma mod- 
eling is implicit pi(J 37 l 38 l ; which numerically damps un- 
resolved, high-frequency modes, and is utilized in codes 
such as LSP^ and ELIXIRiP^. Of course, the reduced- 
model approach is inapplicable to laser-plasma interac- 
tions, or low-density regions with, e.g., Debye sheaths. 
Ion dynamics is not modeled in Zuma (although includ- 
ing them is consistent with the reduced- model approach) , 
and we assume charge neutrality: n e 6 = Zni where n e 6 is 
the number density of free (not atomically bound) back- 
ground electrons and nj = rti is the total ion density. 

Zuma advances each fast electron by 



dx/dt 



dp/dt ■ 



(i + nfi) -v d p + R (5) 
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FIG. 5. Average (solid black) and rms (solid blue, slightly 
higher) polar angle for the angle spectrum in Eq. J3|, vs. the 
parameter AO. The red dashed line is the approximate form 
(0) = 0.691 AO from Eq. @. 



III. ZUMA-HYDRA INTEGRATED MODELING 

Our transport modeling is done with the hybrid-PIC 
code ZumsP^ coupled to the radiation-hydrodynamics 
code HydrsP^. We describe here the Zuma model in some 
detail, and how it is coupled to Hydra. We briefly discuss 
how we run Hydra, and refer the reader to the extensive 
literature on this well-established code. 



A. Zuma hybrid-PIC code 

Zuma is a parallel, hybrid-PIC code that currently sup- 
ports 3D Cartesian and 2D cylindrical RZ geometry. It 
employs an explicit time-stepping approach, treats the 
fast electrons by a standard, relativistic PIC method, and 
models the background plasma as a collisional fluid. The 



where p = m e c^f3 is the relativistic momentum. The 
term —VdP is frictional drag (energy loss), and the 
Langevin term R represents a random rotation of p which 
gives angular scattering. We uss the drag and scat- 
tering formulas of Solodov and BettP^, and Davies et 
a/! 41 * [ We follow the numerical approach of Lemons^, 
by applying drag directly to Ijol and then randomly rotat- 
ing its direction. Manheimer^ presented a similar colli- 
sion algorithm which acts on Cartesian velocity compo- 
nents. Binary-collision algorithms like that of Takizuka 
and Abe^ have advantages like exact conservation rela- 
tions and can be applied to models like ourPS. They gen- 
erally require, however, the drag and scattering to satisfy 
an Einstein relation and thus have the same "Coulomb 
logarithm," which is not the case for the formulas used 
here. An Einstein relation obtains when both processes 
result from many small, uncorrelated kicks to the test- 
particle momentum. Our angular scattering arises from 
such binary collisions, but the energy loss also contains 
a collective (Langmuir- wave emission) part. The energy 
loss is off all background electrons (free and atomi- 
cally bound), and both types of electrons are treated us- 
ing the Solodov- Davies energy-loss formula. This strictly 
applies to free electrons, or to bound electrons in the 
limit where the density effect dominates 47 . Radiative loss 
is neglected, although it becomes important for high-Z 
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ions and high-energy electrons. Specifically, 
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with cj 2 e : 

classical electron radius. This gives rise to a stopping 
power (ds = \v\dt) of 



dE 
ds 



-m e cjf3is d = -47rra e cX-p-L d . (9) 



The Langevin term R is chosen to give the following 
mean-square change in pitch angle iJj : with respect to 



(( A ^) 2 ) _ snort 



As 



72/34 



L e = lnA-(l/2)(l + ln[2 7 + 6]), 



Lj = lnA-(l/2)(l + /3 2 ). 



(10) 



(11) 



(12) 



A = 2\]j e p/h with Aj3 e = eoT e /n e 6e 2 . is the nuclear 
charge, since screening in partially- ionized ions only af- 
fects very small- angle scatters. 

The fast electron current Jf is deposited onto the spa- 
tial grid, and the background current J5 is found from 
Ampere's law: 



Jb 



-J f X V x B. 



(13) 



The magnetic field is advanced by Faraday's law, d t B 
-V x E. 

The electric field is given by the Ohm's law: 



E 

E c 
Enc ■ 



E c + E NC , 



(14) 



V 



J h -e~ l4 f • VT e , 



en eb 



v e b x 5. 



We follow "notation II" of Ref. \M E c and E NC are, 
respectively, collisional and collisionless effects. J5 is 
the background electron current; if ion currents were in- 
cluded, Jfo in Eqs. (13, 14) should be replaced by the 
(total, electron) current. Our Ohm's law neglects terms 
arising from advection ~ • \7v e b, off-diagonal pressure 
terms, and collisions between background and fast elec- 
trons. The background temperature T e (currently the 
same for both electrons and ions) is updated due to col- 
lisional heating, as well as fast-electron frictional energy 



loss Qf r i C (all of which is deposited as heat, not directed 
flow): 



(3/2)(l + l/Z)n eb d t T e =E C -J b + Q 



fric- 



(15) 



We neglect heat flow (e.g. due to gradients) in Zuma, and 
rely on the coupling to Hydra to provide that physics. For 



the collisional transport coefficients and '/? , we use 
the approach of Lee and More^, but with the numerical 
tables of Ref.09]to account for electron-electron collisions 
and background magnetization. We utilize a modified 
version of Desjarlais' extension to Lee and More 51 , and 
use his extension of Thomas-Fermi theory to find the 
ionization state Z. 



B. Hydra for transport simulations 

This section describes how we run Hydra for our cou- 
pled Zuma-Hydra transport simulations. We run in cylin- 
drical RZ geometry on a fixed Eulerian mesh. The radi- 
ation is modeled with implicit Monte-Carlo photonics, 
and tabulated equation of state and opacity data are 
used. Fusion reactions occur in all initially dense (p > 10 
g/cm 3 ) DT zones, with alpha transport and deposition 
done via multi-group diffusion; no neutron deposition 
is done, although this could lower the ignition energy 
slightly. Electron thermal conduction uses the Lee and 
More model with no magnetic fielcpH Although Hydra 
has an MHD package and the option for magnetized ther- 
mal conduction, we currently do not use these features. 



C. Zuma-Hydra coupling 

The coupling of Zuma and Hydra is as follows. Zuma 
models a subset of Hydra's spatial domain, since the fast 
electron source becomes unimportant far enough from 
the source box. This paper reports results in cylindri- 
cal RZ geometry, and 3D Cartesian results have been 
reported iri^ 2 -'. Data transfer between the codes is per- 
formed via files produced by the LLNL Overlink code^, 
which can interpolate between different meshes. 

The two codes are run sequentially for a set of "cou- 
pling steps" that are long compared to a single time step 
of each code. A coupling step from time to to t\ consists 
of: 

1. Plasma conditions (materials, densities, tempera- 
tures) are transferred from Hydra to Zuma. 

2. Zuma runs for several time steps from to to t\. 

3. The net change in each zone's background plasma 
energy and momentum is transferred from Zuma to 
Hydra, as external deposition rates. 

4. Hydra runs for several time steps from time to to 
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Zuma calculates its own ionization state each timestep, 
and does not use Hydra's value. For the results in this 
paper, we ran both codes for 20 ps when the electron 
source was injected and then ran Hydra for 180 ps to 
study the subsequent burn and ignition. Such a run, uti- 
lizing 24 CPUs for Zuma and 48 for Hydra, takes several 
hours of wall time to complete. 3D runs are much more 
demanding, so 2D runs are envisioned for routine design 
studies. 



IV. IGNITION-SCALE MODELING WITH PIC-BASED 
FAST ELECTRON SOURCE 

The next two sections present results of Zuma-Hydra 
modeling of an idealized ignition-scale, cone-guided tar- 
get. This section considers cases with no initial magnetic 
field, and the next studies imposed field schemes to mit- 
igate source divergence. Table [T] summarizes the runs, 
and Fig. 17 contain RZ plots of the ion pressure and fast 
electron current for several runs. 
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FIG. 6. (Color online) Initial target density in g/cm 3 . The 
red line indicates the source cylinder where fast electrons are 
injected. The 8 g/cm 3 carbon cone is colored blue for clarity. 



A. Simulation setup 

We consider a spherical assembly of equimolar DT 
fuel, relevant for high-gain IFE uses. It is depicted 
in Fig. p] The DT mass density is p[g/cm 3 ] = 10 + 
440expP( J R/70/im) 12 ] where R is the distance from 
(r,z) = (0, 117) pm. This gives, for p > 100, an aerial 
density of pR = 3.0 g/cm 2 and mass m = 0.57 mg. With 
the simple burn-up estimate / = pR/(pR + 6) = 1/3, 
we obtain a total fusion yield (neutron and a) of Y = 
fm • 338 MJ/mg = 64 MJ. Igniting such targets at a rate 
of 16 Hz would provide 1 GW of gross fusion power. A 
flat-tip carbon cone is located ~50 pm to the left of the 
dense DT fuel. The cone density of 8 g/cm 3 (2.3x solid) 
was chosen so that, when fully ionized, the total pressure 
is the same in the cone and 10 g/cm 3 DT. All materials 
are initially at a temperature of 100 eV. 

Simulation parameters were as follows. Both codes 
used a uniform mesh with 1 pm cell size. We leave 
the question of beam-plasma micro-instabilities (e.g., 
resistive filamentatiori 3 -^^ or electro-thermaP^) to fu- 
ture work. The Hydra domain extended (in pm) from 
(r, z) = (0, -50) to (200, 250) while Zuma ran on the sub- 
domain (0, —50) to (120,250). The Zuma timestep is set 
by ensuring no fast electrons cross more than one zone 
per step (cAt < max[Ar, Az]) and resolving the elec- 
tron cyclotron frequency. Since Zuma does not support 
light-wave propagation, there is no light-wave Courant 
stability condition. We used At — 0.5 — 1 fs, which gives 
cAt/Az = 0.15 - 0.3 and uj ce At = 1 for B = 57 - 114 
MG (u ce = eB/m e is the non-relativistic cyclotron fre- 
quency). The Hydra timestep was variable, and the cou- 
pling timestep was 0.2 ps for the first 1 ps and 0.5 ps 
subsequently. 

The fast electron source was injected at 



z = —20 pm with an intensity profile If as t(r,t) = 
7 0/ exp [-(l/2)(r/r spot ) 8 ] f(t) with f(t) a flattop from 
0.5 to 18.5 ps with 0.5 ps linear ramps. Unless specified, 
we use r spo t = 18 pm. The total injected fast electron 
energy, £7f ast , is proportional to 2o/ 5 hi particular, 
for our runs, £f ast [kJ] = 7 //5.77 x 10 18 W/cm 2 . As 
discussed in Sec. [TT , we assume Jo/ is 0.52 times the laser 
intensity ioL, which we need to find the ponderomotive 
temperature T p and energy spectrum. We consider a 
527 nm (2nd harmonic of Nd:glass) laser wavelength, 
since this lowers T p oc Ao compared to first harmonic 
light (but is technologically more challenging). 



B. Results with artificially-collimated source 

The fusion yield for the PIC-based energy spectrum 
and an artificially collimated source (AO = 10°) is plotted 
against £f as t for several values of r spot in Fig. [7[ The 
points lie on a somewhat universal curve. This is due to 
two competing effects, both of which are discussed later 
in this section. First, the hot spots are in the "width 
> depth" regime^, where increasing the hot spot radius 
raises the required deposited heat for ignition. On the 
other hand, increasing the source area for fixed power 
decreases the energy of individual electrons, and leads to 
more effective stopping in the hot spot depth. We do not 
expect a strong dependence on r spo t for situations where 
source divergence has been mitigated. We use r spot = 
18 pm in subsequent runs, since this ignites for the lowest 

£ fast Of 132 kJ. 

For a collimated source, an estimate of the minimum 
ignition energy that must be delivered to the hot spot 
is given by Atzeni et al. in Ref. 56. They performed 
2D rad-hydro simulations of idealized, spherical fuel as- 
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Case 


Description 


il/fast lOW 
[KJJ 


i^fast high 
[KJJ 


Yield low 
n\/r ti 

[MJJ 


Yield high 
n\/r ti 

[MJJ 


DQO 


E =1.5 MeV, A0 = 0, no ang. scat, or E/B 


15.8 


18.5 


0.217 


58.1 


DQlO^nono 


E =1.5 MeV, AO = 10°, no E/B 


25.4 


30.4 


0.189 


56.8 


DQ10_noEB 


PIC-based dN/dE, AO = 10°, no E/B 


81.0 


102 


0.928 


54.9 


DQ10 


AO - 10°, B z = 0, full Ohm's law 


121 


132 


0.426 


48.7 


DQ90 


AO = 90°, B z = 


949 


N/A 


6.82E-4 


N/A 


DQ90_36 


DQ90 but r spo t = 36 pm 


1270 


N/A 


0.0144 


N/A 


BZ30 


B z — 30 uniform 


211 


237 


0.538 


52.3 


BZ50 


B z — 50 uniform 


106 


132 


2.66 


54.0 


BZ30-75 


B z = 30 - 75 


316 


N/A 


0.0523 


N/A 


BZ50-75 


£ z = 50 - 75 


158 


211 


0.412 


53.9 


BZO-50 


B z — 0.1 — 50 quickly in z 


211 


N/A 


0.0785 


N/A 


BZ30pipe 


B z = 30 hollow pipe 


290 


316 


0.276 


48.5 


BZ50pipe 


B z = 50 hollow pipe 


132 


158 


0.532 


52.4 


BZ50pipeA 


BZ50pipe but thinner pipe 


185 


211 


0.378 


49.4 



TABLE I. Ignition properties for various Zuma-Hydra cases, all for the initial plasma conditions shown in Fig. [6] Cases starting 
with DQ have no initial magnetic field, and the number indicates the divergence parameter AO in degrees. Those starting with 
BZ use AO — 90° and have an initial axial magnetic field, with the numbers related to the strength in MG. All runs except 
DQ9CL36 have r spo t = 18 pm. The first three runs have no E or B fields, while the others do and use the full Ohm's law. The 
(low, high) £f as t are, respectively, the (largest, smallest) energy which (did not, did) ignite. 
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FIG. 7. (Color online) Fusion yield vs. total injected fast 
electron energy, for Zuma-Hydra runs with an artificially col- 
limated source A# = 10°. r spo t = 10 /xm for black squares 
with solid line, 14 /mi for red circles, 18 /xm for blue triangles 
(case DQ10), and 23 /xm for green crosses. In this and subse- 
quent plots, the dashed black line at 64 MJ is the ideal fusion 
yield described in the text. The blue triangle with .Efast = 132 
kJ is the lowest value that we deem to have ignited. 



semblies heated by a cylindrical beam of mono-energetic, 
forward-going particles which fully stop in a prescribed 
penetration depth. The result is 



£i g = E opt F(pAz,pr), 



opt 



140 



lOOg/cni 



-1.85 



kJ. 



(16) 
(17) 
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FIG. 8. Fusion yield vs. total injected fast electron energy, for 
Zuma-Hydra runs for artificially collimated source A# = 10° 
with r spo t = 18 jLtm (triangles, case DQ10), PIC-based source 
divergence A# = 90° with r spo t = 18 /xm (squares, case DQ90), 
and A# = 90° with r spo t = 36 /im (circles, case DQ9CL36). 



F > 1 if the hot spot transverse radius satisfies pr > 0.6 
g/cm 2 ("width > depth" regime), or the depth satisfies 
pAz > pAz opt = 1.2 g/cm 2 . For our peak density of 450 
g/cm 3 , i? pt = 8.7 kJ, and Atzeni finds an optimal hot- 
spot radius of r opt = 14 pm and pulse length of t opt = 15 
ps. 

The minimum £?f as t which ignited in the DQ10 series 
was 132 kJ, which is 15x E opt . We can understand this 
with the simplified runs listed at the top of Table I. First, 
DQO has a perfectly collimated source (AO = 0), a mono- 
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energetic 1.5 MeV spectrum, no angular scattering, and 
no E or B fields. This ignites for 18.4 kJ, or 2.1E opt . 
This reflects our spot shape and temporal pulse being 
larger than Atzeni's optimal values, and deposition in 
the low-density DT and carbon cone (which Atzeni did 
not include); we also did not optimize the 1.5 MeV source 
energy. The series DQlOjnono uses our small but nonzero 
AO = 10° and includes angular scattering, but no E or 
B fields; we now obtain ignition for 3.5-Eopt- We adopt 
the PIC-based energy spectrum in DQICLnoEB but still 
include no E or B fields. This raises the ignition energy 
by another factor of 3.4, or 12E opt . Finally, turning on E 
and B fields with the full Ohm's law costs another 1.3x, 
bringing us to 15E opt . The role of the energy spectrum 
and Ohm's law for an artificially collimated source is dis- 
cussed in more detail irP^. 

We now estimate the effect of the PIC-based energy 
spectrum on ignition energy. Let E stop be the energy de- 
livered by a perfectly collimated beam with Atzeni's opti- 
mal parameters to a hot spot of depth pAz opt . The total 
fast electron energy E iast = al 0f (with a = 7rr^ pt t opt ) 
is controlled by the fast electron intensity Iof. We con- 
sider Ao = 527 nm and only collisional stopping (no an- 
gular scattering) of the fast electrons. The fraction of 
kinetic energy lost by a fast electron of kinetic energy E 
in the hot spot is well fit by / = min(l, Ex^AE") where 
Edt = 1-3 MeV reflects the stopping in the DT hot spot. 
Integrated over our PIC-based energy spectrum, the ratio 
^stop/^fast is approximately fit by 

iW^fast ~ (l+4/A)sr - 48 los = 1.5-10 19 W/cm 2 . 

(18) 

Figure [9] shows how these formulas apply to our PIC- 
based energy spectrum. 

For Iof ^> Iqs , which our runs satisfy, we obtain 



stop 



Ef a 



r0.52r0.48 
aI 0f 2 0S •> 

£1.92 
^stop 



(otlos) 



0.92 * 



(19) 
(20) 



This is very close to what one finds with a 
ponderomotively-scaled energy spectrum by assuming all 
electrons lose Ed T : 



j?2 

^stop 

aI Ss 



I£ s \ 2 = 13.7 GW 



E 



DT 



<E> 



(21) 

Using the values found above for the bracketed quantities 
and A = 527 nm, we find I* s = 1.6 x 10 19 W/cm 2 . This 
is very close to the fitted value given above. The up- 
shot is that, due to partial stopping of fast electrons, the 
required short-pulse ignitor laser energy oc Ef as t scales 
roughly as square of the hot-spot energy. In addition, 
Ef ast can be decreased by raising Iqs, which would hap- 
pen if the electron stopping power were higher than our 
current value (e.g., due to micro-instabilities 57 or N- 
particle correlated stoppimJ^). 



From Eq. (18), achieving E stop = E opt of 8.7 kJ re- 
quires Ef ast = 5.6E st0 p = 49 kJ. This factor of 5.6x is 1.6 
times larger than the 3.4x we found in going from case 
DQlOjnono to DQ10_noEB, which entailed going from the 
mono-energetic to PIC-based energy spectrum. We con- 
jecture this is because DQlOjnono is already sub-optimal 
enough (ignites for 3.5E op t) that we do not suffer the 
largest possible penalty for using the PIC-based energy 
spectrum. This implies more idealized targets like DQ0 
would pay closer to the full penalty of 5.6x. 



C. Results with PIC-based, divergent source 

Figure [8] presents the results for the cases DQ90 
(^spot = 18 /im) and DQ9CL36 (r spo t = 36 fim) with the 
PIC-based source divergence AO = 90°, as well as case 
DQ10 with an artificially collimated source of AO = 10°. 
The PIC-based source is far from igniting, even for elec- 



tron source energies > 1 MJ. Figure [17| shows the greatly 
increased current divergence for DQ90. Note also the fil- 
aments that develop at large radius. Their nature and 
effect on beam propagation and fuel coupling should be 
further examined in future work. With the realistic di- 
vergence of AO = 90°, the yield for the same Ef as t is 
higher for the larger spot radius. This is generally the 
case for divergent sources, where the benefit of reduced 
laser intensity (and lower energy electrons, which stop 
more efficiently in the hot spot) outweighs the cost of 
increased spot size. Divergent sources thus ignite in the 
so-called "width > depth regime", where the hot spot 
has pr above the optimal value of 0.6 g/cm 2 . The source 
spot size that minimizes Ef as t depends on details, like the 
cone-fuel standoff distance. In any case, its value will be 
unacceptably large for reactor purposes, so we turn our 
attention below to mitigating source divergence. 



V. FAST ELECTRON CONFINEMENT WITH IMPOSED 
MAGNETIC FIELDS 

We now attempt to recover the artificial-collimation 
132 kJ ignition energy, with the realistic source diver- 
gence, by imposing various initial magnetic fields. This 
can be achieved with an axial field B z with no axial vari- 
ation and strength ~50 MG. However, axial variation in 
B z leads to a radial field B r and a v x B force in the 
z direction (i.e., magnetic mirroring), as well as finite 
standoff distance from the source region to the confining 
field. We find that mirroring greatly reduces the bene- 
fit of magnetic fields. A magnetic pipe, with a B z that 
peaks at finite radius, does not suffer from the mirroring 
problem. 

We wish to specify an arbitrary B z (r,z) in cylindri- 
cal coordinates, with no dependence on azimuth <j>. This 
can be accomplished by a vector potential A = A^r, z)(j), 
which by construction satisfies the Coulomb gauge condi- 
tion V • A = 0. The magnetic field B = V x A automati- 
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cally satisfies V B = 0. In particular, B z = (l/^d^rA^) 
and B r = —d z A&. This allows us to solve for As : 



A 4> = - [ dr'r'B z (r',z). 
r Jo 



(22) 



As r 0, ifB z kr p then A^ k(p + 2)~V +1 . Since 
B r scales with r like as long as p > —1, Eq. (22) 



guarantees B r (r — 0) = 0. This is physically necessary, 
since the radial direction is ill-defined at r = 0. We 
can find the current J = J^cj) needed to maintain the 
magnetic field from Ampere's law without displacement 
current: 



Mo^</> = d z B r — d r B z = 



(23) 



We estimate the magnitude of ~ B/noL, and taking 
5-10 MG, L ~ 10 /im gives J - 10 14 A/m 2 . The fast 
electron current is of order which for 527 nm light 

is 1.9 x 10 1T A/m 2 . Although substantial, the currents 
implied by our imposed fields are much less than the 
fast electron current. They may compete with the much 
smaller net (fast plus background) current. 



A. High B z in fast electron source region 

We utilize initial magnetic fields of the form 

B z = B z0 + (B zl - B z0 )H(z) exp [- ((r - r )/Ar) 

(24) 

Except for the magnetic pipe configurations discussed be- 
low, ro = and Ar = 50 /im. We first consider a B z 



<D 
>> 
CZ 

o 
"(/) 

=5 




10 +2 -> 



10 +0 - 



10-* h 



10- 4 - 



Efast = fast e - energy [kJ] 

FIG. 10. (Color online) Fusion yield for runs with r spo t = 
18 /xm. Black solid line with squares is for an artificially col- 
limated source A# = 10° (case DQ10). All other cases use the 
PIC-based source divergence (A0 = 90°). Black dashed line 
with circles is for no imposed B field (case DQ1CL18). The 
other cases have a "uniform" initial B z given by Eq. (24), 
with B z i = 10 MG (red triangles), 30 MG (blue crosses, case 
BZ30), and 50 MG (green diamonds, case BZ50). 



that we call "uniform," since it does not vary between 
the source region and dense fuel. B z (z) for several cases 
is plotted in Fig. 11 For the uniform case, H(z) = 1 
for z < 80 /im, with a piecewise-parabolic ramp to zero 
for z > 110 /im. B z q = 0.1 MG is the uncompressed seed 
field, and we vary the peak compressed field B z \. Setting 
H = 1 slightly decreases the ignition energy, but may be 
less realistic than our H(z). Figure 10 shows the fusion 
yield for the PIC-based source divergence (AO = 90°) 
and various values of B z \. An initial field of B z \ = 10 
MG gives better coupling than the unmagnetized cases, 
but would still lead to an unacceptable ignition energy. 
A field of 30 MG (case BZ30) gives about 2x the ignition 
energy as an artificially collimated source (AO = 10°), 
while 50 MG (case BZ50) gives essentially the same cou- 
pling. 

The pressure and fast electron current profiles in 



Fig. [TT| illustrate the improvement due to the imposed 
field. The current is much more confined in the case BZ50 
than without the field, although not as much as in the 
artificial-collimation case DQ10. The loss of confinement 
at z ~ 100 /im is due to the end of the high-field region 
(see Fig. [TT]). Note also the appearance of current to the 
left of the fast electron injection plane at z = —20 jam. 
We call this the reflected current. It is due to the plasma 
being slightly diamagnetic, and reducing the imposed B z 
somewhat during the course of the run. The reflected 
current becomes enhanced in runs with significant mir- 
roring, such as case BZO-50. 
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We now turn to the effect of more realistic initial field 
geometries. It is plausible to compress the field to the de- 
sired strength, >50 MG, in a fast-ignition fuel- assembly 
implosion 59 . However, it will not be uniform. In particu- 
lar, to the extent the MHD frozen-in law is followed, the 
axial field compression will follow the radial compression 
of matter. Standard schemes of fuel assembly around a 
cone tip will thus result in the largest field being located 
between the cone tip and dense fuel. Moreover, the pur- 
pose of the cone is to provide a plasma-free region so the 
short-pulse laser converts to fast electrons near the fuel. 
The shell motion down the outer cone surface launches a 
strong shock in the cone, which must not reach the inner 
cone surface before the short-pulse laser fires (to avoid 
a rarefaction that would fill the cone interior). In stan- 
dard schemes, the short-pulse laser thus converts to fast 
electrons in a region with essentially the uncompressed, 
seed magnetic field. The field may be enhanced some- 
what by resistive diffusion of compressed field into and 
through the cone material, but we expect the cone and 
surrounding DT to be sufficiently conducting to prevent 
significant diffusion. 

The upshot is the fast electrons must transit from their 
birth region of low field to a region of high field in front of 
the cone. This poses two separate challenges. First, the 
fast electrons may be reflected axially by the magnetic 
mirror effect. Also, they must travel a finite standoff 
distance before the compressed field can impede their 
radial motion. 

We consider the role of mirroring with no standoff, 
by modifying H(z). We first increase B z to 75 MG at 
z = 30 /im (located between the cone tip and the dense 
fuel) while keeping B z fixed at 30 MG (case BZ30-75) or 
50 MG (case BZ50-75) in the source region. The new B z 
profiles are plotted in Fig. [TT] Table I and Fig. [12] show 
that the energy needed to ignite increases slightly for 
case BZ50-75 but substantially for case BZ30-75. This 
demonstrates the significant impact of mirroring for mod- 
est (1.5-2.5x) increases in B z between the source and fuel 
regions. 

To demonstrate further the impact of mirroring, we 
plot in Fig. 13 the reflected fraction, or the ratio of the 



fast electron energy reaching the left edge of the domain 
(z = —50 /im) to Efast- For fixed i£f as t 5 the reflected frac- 
tion is quite small for the uniform field profiles, but sub- 
stantial for the non-uniform ones. The reflected fraction 
actually understates the effect of mirroring. The low- 
energy electrons deposit more of their energy in the hot 
spot, and are also more magnetized (and thus more likely 
to mirror). So, the reflected electrons would have more 
effectively heated the hot spot than typical electrons. 
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FIG. 11. (Color online) B z profiles at r = 0. Thick black 
dash: case BZ30, thick red dash: case BZ50, black solid: case 
BZ30-75, red solid: case BZ50-75, blue solid: case BZO-50. 
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FIG. 12. (Color online) Fusion yield for runs from Fig. 11 
The curves have the same meaning as in that figure, solic 
green X's is case BZ50pipe. The other curves have the same 
meaning as in Fig. [TT] black dash squares: case BZ30, red 
dash crosses: case BZ50, black solid triangles: case BZ30-75, 
red solid circles: case BZ50-75, blue solid diamonds: case 
BZO-50. 



B. Low B z in fast electron source region: the magnetic 
pipe 

We now turn to the situation where the fast electrons 
are born in the uncompressed seed field. First we con- 



sider case BZO-50, where the field rises quickly in z, so 
that standoff is minimized. Figure [TT] shows the B z pro- 
file. The electrons are still subject to the mirror force, 
which results in an ignition energy of £f as t > 211 kJ (blue 
curve in Fig. 12). Runs with higher Ef ast encountered 



numerical difficulties, which we are studying. Figure 13 
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FIG. 13. (Color online) Reflected fraction, i.e. the fast elec- 
tron energy reaching left z boundary, divided by Efast- The 
runs and curves are as in Fig. [12] solid green X's: case 
BZ50pipe, black dash squares: case BZ30, red dash crosses: 
case BZ50, black solid triangles: case BZ30-75, red solid cir- 
cles: case BZ50-75, blue solid diamonds: case BZO-50. 



depicts the reflected fraction, indicating substantial mir- 
roring in this case. Figure [TT| shows the increase in the 
reflected current in case BZO-50 compared to BZ50. 

To remedy mirroring, we propose a hollow magnetic 
pipe, which is free of high field at small radius. Fast 
electrons are reflected by the pipe as they move out- 
ward radially, but do not experience a mirror force in 
z (d z B z — B r = inside the pipe). A certain product 
of field strength times length is needed to reflect an elec- 
tron, and can be estimated for planar (not cylindrical) 
geometry by^ 



BL > K 7 /3(l-cos<9), 



K 



17.0MG-/im. 

(25) 

For our PIC-based AO — 90° angle spectrum, 
(1 — cos#) = 0.43. The ignition energy for case DQ10 
(artificial collimation) was 132 kJ, which gives an aver- 
age electron energy of 8.5 MeV. This requires BL > 129 
MG-/im to reflect. Although lower-energy particles are 
easier to reflect and stop more fully in the DT hot spot, 
our spectrum does not contain much energy there. For 
B = 50 MG, we need a pipe of thickness 2.6 jam. We con- 
sider pipes that are thicker than this, which substantially 
reduce the ignition energy over the no-field case. The re- 
sults shown here establish the feasibility of the pipe. We 
are exploring thinner, more optimal pipes, which entails 
variation of other parameters like spot size. 

The initial field for the case BZ50pipe is given by 
Eq. (24) with ro = 35 /mi and Ar = 15 /im, and is dis- 
played in Fig. 14 This gives roughly the same ignition 
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FIG. 14. (Color online) Initial B z profile, in MG, for case 
BZ50pipe. 



energy as the uniform B z field, shown as the solid green 



line in Fig. 12 The mirror effect, as measured by the 
reflected fraction, is about the same as for the uniform 
case BZ50 (see Fig. [l3]). We studied pipes with a peak 
field of 30 and 50 MG (cases BZ30pipe and BZ50pipe), 
as well as thin and thick 50 MG pipes (cases BZ50pipe 
and BZ50pipeA, see Fig.[l6|). 

We summarize the development of this paper in 
Fig. 15 The challenge was to find imposed magnetic 
field configurations that recover the performance of an 
artificially collimated fast electron source, when using the 
PIC-based divergent source. A uniform 50 MG axial field 
does this, and may even perform slightly better. How- 
ever, the more realistic case is for fast electrons to be 
born in a lower field and suffer magnetic mirror forces. 
To circumvent this, we introduced the hollow magnetic 
pipe. For a 50 MG peak field, this works essentially as 
well as the uniform field. A lower peak field of 30 MG 
performs significantly worse than the 50 MG cases, for 
both the uniform and pipe configurations. 



VI. CONCLUSION 

In this paper, we have presented recent transport mod- 
eling efforts geared towards a fast ignition "point design." 
This requires knowledge of the fast electron source pro- 
duced by a short-pulse laser. We characterized the results 
of a full-PIC 3D simulation with the PSC code in terms 
of ID energy and angle spectra. The energy spectrum is 
well-matched by a quasi two-temperature form, which we 
scale ponderomotively as we vary IX 2 . The angle spec- 
trum is divergent, and the PIC data showed only a slight 
reduction in average angle with electron energy (justify- 
ing our ID factorization). A more sophisticated handoff 
method involving a 4D distribution function, and using 
more recent PIC simulations, seems to give qualitatively 
similar resultPl The major design challenges posed by 
this source are: 1. the electrons are too energetic to fully 
stop in a DT hot spot, and 2. they are sufficiently di- 
vergent that mitigation strategies are required in a point 
design. 
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FIG. 15. (Color online) Fusion yield for run cases DQ10 (solid 
black squares), BZ50 (solid red triangles), BZ50pipe (dashed 
red diamonds), BZ50pipeA (dashed blue X's) BZ30 (solid green 
crosses), and BZ30pipe (dashed green circles). 
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FIG. 16. (Color online) Radial B z profiles for pipe cases 
BZ30pipe and BZ50pipe (solid black), and BZ50pipeA (thick 
dashed blue). 

We have developed a transport modeling capability 
which entails the hybrid-PIC code Zuma and rad-hydro 
code Hydra running in tandem. We detailed the physics 
contained in Zuma. It is similar to other codes that 
use a reduced model to eliminate light and Langmuir 
waves. Namely, the displacement current is removed from 
Ampere's law, and the electric field is found from Ohm's 
law (obtained from the background electron momentum 
equation). This model is applicable in sufficiently colli- 
sional plasmas, and for time and space scales longer than 
the plasma frequency and Debye length. 

Zuma-Hydra 2D cylindrical RZ runs on an idealized 
cone-fuel assembly were performed. For a perfectly par- 
allel source (AO = 0, no angular scattering) of mono- 



energetic 1.5 MeV electrons, ignition occurred for 18.5 
kJ of fast electrons, or 2.1x Atzeni's ideal estimate of 8.7 
kJ. We discussed the impact of (small) divergence, E and 
B fields, and the PIC-based energy spectrum, and refer 
the reader to Ref. |22]for more details, including the role 
of different terms in Ohm's law. With the PIC-based 
energy spectrum, the full Ohm's law, and an artificially- 
collimated source (AO = 10°), the ignition energy was 
raised to 132 kJ, or 15x the ideal value. 

The realistic angular spectrum was then considered, 
and raised the ignition energy to > 1 MJ. Several mitiga- 
tion ideas have been proposed, including magnetic fields 
produced by resistivity gradients (e.g. at material inter- 
faces). While this approach is promising, we chose to 
examine imposed axial magnetic fields. An initial, uni- 
form field of 50 MG recovered the 132 kJ ignition energy 
of the artificially collimated source. Assembling ~ 10's 
MG field strengths in an ICF implosion, via the frozen-in 
law of MHD, is reasonable, and has been demonstrated 
recently at Omega. 

However, a cone-in-shell implosion is not likely to pro- 
duce a uniform magnetic field. In particular, the field in 
the fast electron source region (inside the cone tip) will 
not be enhanced much over the seed value, but would be 
enhanced in the region between the cone tip and dense 
fueP^. Fast electrons would therefore encounter an in- 
creasing axial field, and be subject to magnetic mirror- 
ing. Simulations that quantify this effect for a few profiles 
were shown. We showed one way to provide confinement 
but avoid mirroring is a magnetic pipe, which peaks at a 
finite radius. 

We have started to address the design problem of as- 
sembling a pipe field in an implosion. Inserting an ax- 
ial structure (a "wire") between the cone tip and fuel, 
that does not get compressed, is one way to achieve this. 
Magnetic confinement schemes based on self-generated 
azimuthal fields due to resistivity gradients require a 
similar structure. Thus both approaches share some hy- 
dro assembly features, and would mutually benefit from 
progress in hydro design. One advantage of the pipe in 
this regard is that the resistivity (e.g., Z) is irrelevant 
for the pipe, as long as it isn't compressed and thereby 
produce a large on-axis field. The self-azimuthal fields 
rely on resistivity gradients, usually achieved by a high- 
Z material on-axis. This can lead to unacceptable fast- 
electron energy loss or angular scattering in the wire^M 
Megagauss magnetic fields may lower the ignition thresh- 
old by reducing electron thermal conduction out of the 
hot spot or, at even higher values, enhancing alpha de- 
position. 

Integrated hybrid-PIC and rad-hydro simulations offer 
a powerful new tool for fast-ignition modeling, and we 
look forward to them enabling the emergence of attrac- 
tive ignition designs. 
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FIG. 17. (Color online) Ion pressure (left) and fast electron current density | J f as t| (right) at time 10 ps (middle of fast electron 
time pulse) for cases, from top to bottom: DQ10, DQ90, BZ50, BZO-50, and BZ50pipe. All cases have E{ ast = 158 kJ, except 
DQ90 has E iast = 317 kJ. 
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